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ABSTRACT. The first prime number with the special property that 
its addition with its reversal gives as result a prime number too is 
299. The prime numbers with this property will be called Luhn 
prime numbers. In this article we intend to present a performing 
algorithm for determining the Luhn prime numbers. Using the pre- 
sented algorithm all the 50598 Luhn prime numbers have been, for 
p prime smaller than 2 • 10 7 . 



1. Introduction 

The number 299 is the smallest prime number that added with his 
reverse gives as result a prime number, too. As 1151 = 229 + 922 is 
prime. 

The first that noted this special property the number 229 has, was 
Norman Luhn (after 9 February 1999), on the Prime Curios website 
|j2fl . The prime numbers with this property will be later called Luhn 
prime numbers. 

In the Whats Special About This Number? list D3J, a list that con- 
tains all the numbers between 1 and 9999; beside the number 229 is 
mentioned that his most important property is that, adding with his 
reversal the resulting number is prime too. 

The On-Line Encyclopedia of Integer Sequences, [6} A061783], presents 
a list 1000 Luhn prime numbers. We owe this list to Harry J. Smith, 
since 28 July 2009. On the same website it is mentioned that Harvey 
P. Dale published on 27 November 2010 a list that contains 3000 Luhn 
prime numbers and Bruno Berselli published on 5 August 2013 a list 
that contains 2400 Luhn prime numbers. 

2. Smarandache's function 

The function p : N* — * N*, pin) = m, where rri is the smallest 
natural number with the property that n \ m\ (or m\ is a multiple 
of n) is know in the specialty literature as Smarandache's function, 
flZH8]]. The values resulting from n = 1, 2, . . . , 18 are: 1, 2, 3, 4, 5, 3, 
7, 4, 6, 5, 11, 4, 13, 7, 5, 6, 17, 6. These values were obtained with 
an algorithm that results from p's definition. The program using 
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this algorithm cannot be used for n > 19 because the numbers 19!, 
20!, . . . are numbers which exceed the 17 decimal digits limit and the 
classic computing model (without the arbitrary precisions arithmetic 
IfTOlP will generate errors due to the way numbers are represented in 
the computers memory 

3. KEMPNER'S ALGORITHM 

Kempner created an algorithm to calculate p(n) using classical fac- 
torization n = p v f ■ 2 ■■■ p% 3 , prime number and the generalized nu- 

meration base (oj)[ Pi ], for i = l,s, 0]. Partial solutions to the al- 
gorithm for p(n)'s calculation have been given earlier by Lucas and 
Neuberg, 0. 

Remark 3.1. If n e N*, n can be decomposed in a product of prime 
numbers n = pf 1 ■ p ^ 2 ■ ■ • p" s , were p, are prime numbers so that pi < 
P 2 < . . . < p s , and s > 1, thus Kempner's algorithm for calculating 
the /i function is. 



Mn.) = max | Pl ■ ) (pi) . V- ' . 

P *' ’ 

where by (^[p]) (pj we understand that a is " written " in the numer- 
ation base p (noted n^) and it is "read” in the p numeration base 
(noted /3(p), were (3 = a^), JH p. 39]. 

4. Programs 



The list of prime numbers was generated by a program that uses 
the Sieve of Eratosthenes the linear version of Pritchard, H5J, which 
is the fastest algorithm to generate prime numbers until the limit of 
L, if L < 10 s . The list of prime numbers until to 2 • 10' is generated 
in about 5 seconds. For the limit L > 10 8 the fastest algorithm for 
generating the prime numbers is the Sieve of Atkin, |[T|. 



Program 4.1. The Program for the Sieve of Eratosthenes, the linear 
version of Pritchard using minimal memory space is: 



CEPbm(L) := 



A A- floor (|) 
for k e 1..A 
is-primek 1 
prime <(—(2 3 5 7) T 
i <(— last (prime ) + 1 
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for J e 4, 7.. A 
is -prime j 0 

k <- 3 

s 4— (prime k~ i) 2 
t ( prime k ) 2 

while t < L 

for j Et,t + 2 • primek-.L 

is-primej-i <— 0 

for j E s + 2, s + 4.T — 2 

if is -prime j-i = 1 
2 

primes j 
i •*— i + 1 

•S 4 — t 

k <— k + 1 

t (prime k) 2 

for j G s + 2, s + 4..L 

if is jprime j -i=l 
2 

primei 4— j 
i •<— i + 1 
return prime 

Program 4.2. The factorization program of a natural number; this 
program uses the vector p representing prime numbers, generated 
with the Sieve of Eratosthenes. The Sieve of Eratosthenes is called 
upon trough the following sequence: 

L 2 ■ 10' t 0 = time( 0) p := CEPbm(L) t\ = time( 1) 

(fi — f 0 ) s = 5.064s last(p) = 1270607 pi as t { P ) = 19999999 

Fa(m ) := return (”m = ” m ” > ca ultimul p 2 ”) if m > ( Pi as t(p )) 2 
J <-l 
k <- 0 

while m > p^ 

if mod (m,pj )= 0 

k i — fc 4 1 
m 

mi 

Pj 

otherwise 

f <— stack[f , (pj, k)] if k > 0 
J <- J + 1 
k <- 0 
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/ A- stack[f, ( pj , k)] if k > 0 
return submatrix(f , 2, rows(f), 1, 2) 

We presented the Kempner's algorithm using Mathcad programs 
required for the algorithm. 

Program 4.3. The function counting all the digits in the p base of nu- 
meration in which is n. 

ncb(n,p ) := return ceil(\og(n,p)) if n > 1 
return 1 otherwise 

Where the ceil(-) Mathcad function represents the upper non-decimal 
number. 

Program 4.4. The program intended to generate the p generalized 
base of numeration (noted by Smarandache [p]) for a number with 
m digits. 

a(p,m) := fori G l..m 

pi _ i 

(k < — 

p — 1 

return a 

Program 4.5. The program intended to generate for the p base of nu- 
meration (noted by Smarandache (p)) to write the a number. 

b(a,p ) := return (1) if p = 1 
for i G 1 ..ncb(a,p) 
b i k- p 1-1 
return b 

Program 4.6. Program that determines the digits of the generalized 
base of numeration [p] for the number n. 

Nbg(n,p) := m 4— ncb(n,p) 
a <— a(p, m) 
return (1) if m = 0 
for i Em.. 1 

( n 

Ci <— trunc — 

\ a * 

n G- mod (n, af) 

| return c 

Program 4.7. Program for Smarandache's function. 

p,(n) := return ” Err. n nu este intreg ” if n ^ trunc (n ) 
return "Err. n < 1” if n < 1 
return (1) if n= 1 
/ <— Fa(n) 
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p /W 
a <- /< 2 > 

for k = l..rows(p) 

V k^Pk- Nbg(a k ,p k ) ■ b{a k ,p k ) 
return max(^) 

This program calls the Fain) factorization with prime numbers. 
The program uses the Smarandache's [34] Remark - about the Kemp- 
ner algorithm. The p.prn file generation is done once. The reading of 
this generated file in Mathcad's documents results in a great time- 
save. 



Program 4.8. Program with which the file p.prn is generated 



VFp(N) : = 



/ii 1 

for n G 2..N 
p n <— p{n) 
return p 



This program calls the |4~7| program for calculating the value of the p 
function. The sequence of the p.prn file generation is: 

t 0 := time iO) W RIT EPRNi” p.prn'’’) := VFp{ 2 • 10 7 ) t\ := time{ 1) 

(ti — to)sec = ”5 : 17 : 32.625” hhmmss 

Smarandache's function is important because it characterizes prime 
numbers - through the following fundamental property: 

Teorema 4.9. Let be p an integer > A, than p is prime number if and only 
tf Pip) = V- 

Proof. See 0 p. 31]. □ 

Hence, the fixed points of this function are prime numbers (to 
which is added 4). Due to this property the function is used as pri- 
mality test. 

Program 4.10. Program for testing p's primality. This program re- 
turns the 0 value if the number is not prime number and the 1 value 
if the number is a prime. The file p.prn will be read and it will be 
assigned to the p vector. 



ORIGIN := 1 p:= READPRN (” . . . \p.prn ”) 

Tppin) := return ” Err. n < 1 sau n f Z” if n < 1 V n truncfn) 
if n > 4 

I return 0 if p n ^ n 
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| return 1 otherwise 
otherwise 

return 0 if n= 1 V n= 4 
return 1 otherwise 

Program 4.11. Program that provides the reveres of the given m num- 
ber. 

R(m ) := n floor(\og(m)) 
x A- rn ■ 10 _n 
for k G 1 ..n 
Ck <— truncfx ) 
x <— (x — Cfc) • 10 
c n+ i roundfx ) 

Rm <— 0 
for k E n + 1..2 
-Rm (Rm + Cfc) • 10 
return Rm + ci 

Program 4.12. Search program for the Lw/m prime numbers. 

PLuhn(L) := n A- last(p ) 

5 G- (229) 

A; •< — 51 
while pk < L 
N 4— R(p k ) +p k 
S <— stack(S,pk ) if Tpp(N) = 1 
k <— k + 1 
return S 

The initialization of the S stack is done with the vector that con- 
tains the number 229. The variable k has the initial value of 51 be- 
cause the index of the 229 prime number is 50, so that the search for 
the Luhn prime numbers will begin with p 5i = 233. 

5. List of prime numbers Luhn 

We present a partial list of the 50598 Luhn prime numbers smaller 
than 2 ■ 10 7 (the first 321 and the last 120): 

229 239 241 257 269 271 277 281 439 443 463 467 479 499 613 641 653 
661 673 677 683 691 811 823 839 863 881 20011 20029 20047 20051 
20101 20161 20201 20249 20269 20347 20389 20399 20441 20477 20479 
20507 20521 20611 20627 20717 20759 20809 20879 20887 20897 20981 
21001 21019 21089 21157 21169 21211 21377 21379 21419 21467 21491 
21521 21529 21559 21569 21577 21601 21611 21617 21647 21661 21701 
21727 21751 21767 21817 21841 21851 21859 21881 21961 21991 22027 
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22031 22039 22079 22091 22147 22159 22171 22229 22247 22291 22367 
22369 22397 22409 22469 22481 22501 22511 22549 22567 22571 22637 
22651 22669 22699 22717 22739 22741 22807 22859 22871 22877 22961 
23017 23021 23029 23081 23087 23099 23131 23189 23197 23279 23357 
23369 23417 23447 23459 23497 23509 23539 23549 23557 23561 23627 
23689 23747 23761 23831 23857 23879 23899 23971 24007 24019 24071 
24077 24091 24121 24151 24179 24181 24229 24359 24379 24407 24419 
24439 24481 24499 24517 24547 24551 24631 24799 24821 24847 24851 
24889 24979 24989 25031 25057 25097 25111 25117 25121 25169 25171 
25189 25219 25261 25339 25349 25367 25409 25439 25469 25471 25537 
25541 25621 25639 25741 25799 25801 25819 25841 25847 25931 25939 
25951 25969 26021 26107 26111 26119 26161 26189 26209 26249 26251 
26339 26357 26417 26459 26479 26489 26591 26627 26681 26701 26717 
26731 26801 26849 26921 26959 26981 27011 27059 27061 27077 27109 
27179 27239 27241 27271 27277 27281 27329 27407 27409 27431 27449 
27457 27479 27481 27509 27581 27617 27691 27779 27791 27809 27817 
27827 27901 27919 28001 28019 28027 28031 28051 28111 28229 28307 
28309 28319 28409 28439 28447 28571 28597 28607 28661 28697 28711 
28751 28759 28807 28817 28879 28901 28909 28921 28949 28961 28979 
29009 29017 29021 29027 29101 29129 29131 29137 29167 29191 29221 
29251 29327 29389 29411 29429 29437 29501 29587 29629 29671 29741 
29759 29819 29867 29989 . . . 

8990143 8990209 8990353 8990441 8990563 8990791 8990843 8990881 
8990929 8990981 8991163 8991223 8991371 8991379 8991431 8991529 
8991553 8991613 8991743 8991989 8992069 8992091 8992121 8992153 
8992189 8992199 8992229 8992259 8992283 8992483 8992493 8992549 
8992561 8992631 8992861 8992993 8993071 8993249 8993363 8993401 
8993419 8993443 8993489 8993563 8993723 8993749 8993773 8993861 
8993921 8993951 8994091 8994109 8994121 8994169 8994299 8994463 
8994473 8994563 8994613 8994721 8994731 8994859 8994871 8994943 
8995003 8995069 8995111 8995451 8995513 8995751 8995841 8995939 
8996041 8996131 8996401 8996521 8996543 8996651 8996681 8996759 
8996831 8996833 8996843 8996863 8996903 8997059 8997083 8997101 
8997463 8997529 8997553 8997671 8997701 8997871 8997889 8997931 
8997943 8997979 8998159 8998261 8998333 8998373 8998411 8998643 
8998709 8998813 8998919 8999099 8999161 8999183 8999219 8999311 
8999323 8999339 8999383 8999651 8999671 8999761 8999899 8999981 
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6. Conclusions 

The list of all Luhn prime numbers, that totalized 50598 numbers, 
was determined within a time span of 54 seconds, on an Intel pro- 
cessor of 2.20 GHz. 
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